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THE COMPRESSIBLE LAMINAR BOUNDARY LAYER WITH HEAT 
TRANSFER AND ARBITRARY PRESSURE GRADIENT 
By Clarence B. Cohen and Eli Reshotko 


SUMMARY 

An approximate method for the calciLLation of the compressible lam- 
inar boimdary layer with heat transfer and arbitrary pressure gradient, 
based on Thraites ' correlation concept, is presented. The method re- 
sults from the application of Stewartson's transformation to Praudtl's 
equations, which yields a nonlinear set of two first-order differential 
equations. These equations are then expressed in terms of dimensionless 
parameters related to the wall shear, the surface heat transfer, and the 
H transformed free-stream velocity. Thwaites ' concept of the unique Inter- 

g dependence of these parameters is assumed. The evaluation of these 

quantities is then carried out by utilizing exact solutions recently 
obtained . 

With the resulting relations, methods are derived for the calcula- 
tion of the two-dimensional and axially symmetric laminar boixndary layer 
with an arbitrary free-stream velocity distribution, Mach number, and 
surface temperatui’e . 

The combined effect of heat transfer and pressure gradient is dem- 
onstrated by applying the method to calculate the characteristics of 
the boundary layer on thin supersonic surfaces and in a highly cooled, 
convergent -divergent, axially symmetric rocket nozzle. 


INTRODUCTION^ 

In recent years ^ with the advent of laminar airfoils and with the 
observation of laminar boundary layers at Reynolds numbers as high as 
50^^10^ (ref. 2 ) ^ the ability to reliably estimate viscous flow and heat- 
transfer effects for a laminar boundary layer has become increasingly 


The principal developments of this paper^ which is part of the 
doctoral dissertation of the senior author (ref. l)^ were carried out 
under the stimulus and guidance of Professor Luigi Crocco and the spon- 
sorship of the Daniel and Florence Guggenheim Foundation. The final 
analysis and the computations were completed at the NACA Lewis laboratory. 
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important. Moreover^ with high-altitude flight becoming more common, 
the subsequent lower Reynolds numbers encountered should more frequently 
produce a laminar boundary layer. Stability calculations based on the 
theory of Lees and Lin (ref. 3) have also emphasized the possibilities 
of maintaining a laminar boundary layer through cooling of aerodynamic 
siirfaces. The effect of favorable pressure gradients in increasing the 
stability of laminar boundary layers may also make solutions to the lam- 
inar problem applicable to the design of nozzles and turbine blades . 

Solution of the laminan -boundary-layer equations which Include ef- 
fects of compressibility, pressure gradient, and heat transfer have been 
quite limited in number. Of the exact solutions, most have restrictions 
of range or application, or both. The solutions of references 4 and 5 
are restricted to zero pressure gradient, while those of reference 6 
allow sma.1 1 pressure gradients. The developments of reference 7 are 
restricted to small heat transfer and low Mach number. Solutions ob- 
tained by assuming that fluid properties are constant or that the Mach 
number is essentially zero are obtained in references 8 to 10. Those 
solutions of references 11 and 12 that are for a Prandtl number of unity 
are not restricted in range of compressibility, pressure gradient, or 
heat transfer. However, they apply to specific types of free-stream 
velocity distribution which are inappropriate for general practical 
problems . 

In 1921, von Karraan (ref. 13) recognized that to solve the skin- 
friction problem it was not necessary to have the exact and complicated 
solution, but that it would be quite satisfactory to evaluate average 
quantities across the layer if they could be related to the surface 
values . The concepts of displacement and momentiam thicknesses were in- 
troduced, thxis considerably simplifying the mathematics of the problem. 
With this integral method, if the form of the velocity profile is re- 
lated to a single parameter, a method of calculating the boundary layer 
is obtained. Pohlhausen (ref. 14) carried out this method by postulat- 
ing a quartic velocity profile depending upon the local pressure gra- 
dient. A number of investigators have extended Pohlhausen 's method to 
compressible flows over insulated surfaces. 

With the presence of heat transfer at the surface, the compressible 
problem becomes more complex. Kalikhman (ref. 15) defined certain heat- 
flow quantities analogoiis to the displacement and momentum thicknesses 
and, in a manner similar to Pohlhausen 's, developed a complex iterative 
procedure for the solution of the general problem. More recently, ref- 
erences 16 to 19 have further developed this technique. The preceding 
methods are tedious, since they require a solution of at least one or- 
dinary differential equation for any particular problem. 

Thwaites' method (ref. 20 ) does not require the solution of ordi- 
nary differential equations. In that formulation, it is suggested that 
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the basic goal of an integral approach might be achieved if the nrohlPTn 
IS considered as that of relating the shear, its nor^ dLS^ve 
at the vallj, and the form factor (ratio of displacement thickness) to 
one another without specifying a type of profile. To thifeS nLS- 
mensional fonts of these quantities were defined and were evSlatS W 

tl'e incompressible laminar boundaiy layer 
auanStler-f nearly universal relation existed between tSse 

quantities for favorable pressure gradients and, for adverse pressure 
gradients, Thwaites selected a single representative relation!^ A unique 
correlation was chosen which reduced the solution of an incompressible 
problem to the evaluation of a single integral.^ 

trantthP^f,^^ recognized that, in the absence of heat 

b Vi ^ I’^^ndtl number of unity, the Ulingworth-Stewartson 
correlation between compressible and incompressible boundary -layer solu- 
tions Crof. 23) could be used to extend Thwaites' method tfiiSSe e)- 
lects of compressibility. 

With the presence of heat transfer the application of Stewartson's 
transformation does not correlate a given con^ressible problem to an 
equivalent incompressible one. Thus, the universal relation previously 
described IS not adequate. Unfortunately, there is little possibility 
of establishing a family of "universal" relations with, for example, the 
wall ten5>erature as the distinguishing parameter, since a variety of 
exact solutions to this problem are not available. However, one such 
set of relations may be obtained from the solutions of references 11 and 

-L^ m 


In the present paper, after formulation of a nonlinear system of two 
first-order differential equations (with the major restriction being a 
3^ear viscosity law), methods of solution are developed depending on 
Thwaites concept of universal functions. The functions used for this 
puiq)ose are evaluated from the solutions of reference 12 only. 


BOUNDARY-LAYER EQUATIONS 

The equations of the steady, two-dimensional compressible l Rmirm-r 
boundary layer for perfect fluids are 


Other approaches, such as that of Young and Winterbottom (ref. 21 ) 
r expressions for the momentimi thickness similar to that^^ 

ot thwaites. In that analysis, however, the derivation was a modiflca- 
tion of the PohJiiausen technique. The application of a correlation con- 
cept was Aot proposed. 
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Cont inuity: 


Momentiim: 


(pu) + ^ (pv) = 0 


Su , c3u ^ I ^ ^ 


1 = ^ 



(1) 




Energy: 


Sh 


Sh 


u 


"§y \Pr 



AH symbols are defined in appendix A. 
The viscosity law to be assumed is 

ji- = 


P-O 


"^0 


(3) 


(4) 


Equation (i) is of the form taken by Chapman and Rubesin (ref. s), except 
that the reference conditions (Po^'^o) free-stream stagnation values, 

since, ia the case of pressure gradient, the local "external" values are 
not constant along the flow. The constant X is used to imtch the vis- 
cosity with the Sutherland value at a desired location. If this loca- 
tion is assumed to be the surface, the result is 


X = 


^ tp + ^ su 

t^ + kgy_ 




(5) 


where 


k = Sutherland's constant = 216° R (for air) 
su 

Stewartson's transformation. - By using the following definitions 


of the stream function: 


, _ -HH 

V PQ 


( 6 ) 


'If. 


= 

“ pq 
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a slight modification of Stewartson's transformation (ref. 23) mav be 
written 

^■e Pe ^ 

dX = X -H (ix 

pq 

a r 

dY = dy 

PQ % J 

The transformed coordinates are now represented by upper-case letters 
(X^ Y) ^ and the subscript e refers to local conditions at the outer 
edge of the boundary layer (external) where the flow is assumed to be 
isentropic. The subscript 0 refers to free-stream stagnation values. 

Applying equations (l) to (?) to the boundary -layer equations (l) 
to ( 3 ) and assuming that Pr and c^ are constant (but not yet requir- 
ing that Pr = 1 ) result in the following system: 


+ Vy ^ 0 

+ ™Y “ + S) + V(,U^ 


USx + VSx = »o< 


Prf 


where the enthalpy term S is defined for convenience as 


where h^ is the local stagnation enthalpy^ and where the stream func- 
tion has been replaced by the transformed velocities (U, V) through the 
relations 


U = ^Ity 


V = 


The resulting relation between the transformed and physical longitudinal 

aQ 

velocities is U =» — u. 
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The boundary conditions applicable to the system (8)^ (9)^ and (lO) 
for a specified wall temperature are 


u(x,o) = 0 
v(x,o) = 0 
s(x,o) = s^(x) 

lim S = 0 

Y OO 

lim U = U„(X) 

Y ® 



(13) 


Integral equations. - An alternate form of the momentum equation 
may be obtained by subtracting the momentum eqmtion (9) from the prod- 
uct of the continuity equation (8) and the quantity (Ug - U). This 
results in 



XUg - U) + UgUg S = - V 


If equation (l4) is integrated with respect to Y between the limits 
Y = 0 and Y = A, where A is a constant distance normal to the sur- 
face sufficiently large that the conditions S = 0 and U — Ug can 
both be satisfied; there results 

^ (e.u2) + UgUg^(6* + ^) = (15) 


where the transformed momentum thickness 9 ^, the transformed displacement 
■thickness ^i^ and the enthalpy thickness ^ ; are defined as 
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By carrying out the indicated differentiation, equation (l5) 
xn the form ^ ^ 


can be put 




{ 29 ^ + + 






(17) 


eqiiation has the form of the conventional Karman momentum 
with the exception of the enthalpy term 


integral 


_ It should be noted that because of Stewartson’s transformation a 
simple relation exists between the parameter 0. and the actual physi- 
cal momentum thickness 6 . This relation is 


r+1 



(18) 


Following a procedure with the energy equation similar to that for 
the momentiom equation results in 


Ug 

^ ^ F 


_!2_ 

PrU^ \^l 


(19) 


where the enthalpy -flux thickness is defined by 



( 20 ) 


The method presented in this report uses exact solutions to the 
boundary-layer equations including the energy equation. Since both 
the skin-friction and heat -transfer parameters from the exact solutions 
are correlated with a parameter which is evaluated from only the momen- 
tum integral equation^ it will not be necessary to use equations (l9) 
and (20). ^ ^ 


EEDUCED BOUWDARY-IAYER EQUATIONS 

At this stage^ the relation governing the boundary-layer develop- 
ment is equation (l7)^ subject to the boundary conditions ( 9 .) =0 

or where the subscript sp indicates stagnation-point values. 

The former condition on applies when the boundary layer is initiated 
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without a stagnation point (such as in the case of a supersonic tMn 
airfoil). The value of (0^)^^ depends on the value of (Ug^)gp ^ 

on the surface temperature. Values of 0gp are presented in 
talole I- 

Before consideration of a solution which depends on a correlation 
similar to that of Thwaites, it is expedient to transform the preceding 
system of equations to a system involving dimensionless parameters. Tne 
correlation concept will then be introduced and methods of solution 
developed. 


Parametric Form 

The dimensionless parameters, which are related to terms appeartog 
in equations (l7) and (l9), can be defined and evaluated from the follow- 

ing expressions: 


Shear parameter, or first-derivative parameter, 

\^J ~ u t \SyJ 




( 21 ) 


Correlation number (related to pressure gradient), or second-derivative 
parameter, 


U. 


0f VdY^j 


w 


Up G" 


w 


fe)l) 


( 22 ) 


Heat-transfer parameter, or third-derivative parameter, 

w 




(23) 


In the definitions (22) and (23), use is made of the following relations, 
respectively: 




U^U, 


e'^e. 


0 


(1 + 


(24) 


which is obtained from equation (9), and 
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Ue^larV 

A ^ 


(25) 


which is obtained by differentiating equation (9) with respect to Y 
and evaluating the resulting expression at Y = 0. 

20iUe 

If equation (l7) is multiplied by ^ there results 


CVJ ' 

I 

o 


Ue d(0?) 
Vq dX 


= - 2 


, J 0.^= ^ (i)„ 


(26) 


^i + 

where is the form factor for low-speed flow * 0 j . 

Now^ if the definitions (2l) and (22) are inserted^ equation (26) 
becomes 


-U, 


e dX 


fv] “ 2 [“(“me + 2) + j] 


(27) 


A similar procedure can be carried out with the energy equation^ 
although it is not necessary for the calculation method herein presented. 


Correlation in Terms of n 

If^ in the definition of (^1-(26))^ is set equal to zero^ 

equation (26) becomes Thwaites ' momentum equation. If^ in this case_, 
physically valid relations H(n) and l(n) can be established^ the equa- 
tion can be integrated and the problem is solved. The assumption of the 
form of the velocity profile serves this puip)Ose^ and the resultant pro- 
cedure is the well-known Pohlhausen technique. (Kalikhman (ref. 15) 
was the first to carry the same approach over to the case of the thermal 
profile.) Thwaites used the more direct approach by determining whether 
universal relations H(n) and Z(n) could possibly be established from 
the well-known "exact” solutions of the boundary -layer problem. An ex- 
amination of these solutions proved that for favorable pressure gradients 
a single relation for each of these quantities could be established with 
a fair degree of accuracy^ but for adverse pressure gradients the rela- 
tions departed from each other considerably^ as indicated for Z(n) in 
sketch (a) (taken from ref. 20) . Since the correlation technique requires 
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Incompressible flow, insulated surface 



(a) 


the selection of a single relation between all boundary-layer quantities 
regardless of the history of the development of the specific boundary 
layer under consideration, the assiamption of a single relation for all 
boundary layers is not exactly valid. Thwaites chose a relation for 
l{n) that would match the Howarth solution (ref. 24) at separation. The 
Pohlhausen solution predicts separation to occur much later than any of 
the other solutions in the sketch, while the Hartree solution (ref. 25) 
predicts separation to occur earlier than the Howarth solution. Stew- 
antson (ref. 23) has indicated that under certain conditions Howarth 's 
solution would predict separation too late. 

With the exception of the described differences in the relation 
between boundary-layer quantities (due to the selection of the solutions 
of ref. 12 for the evaluation of the correlated quantities), the method 
of correlation to be presented herein will contain as special cases the 
method of Thwaites for incompressible flow and the results of Rott and 
Crabtree for compressible flow over Insulated surfaces. 

The concept of correlation is herein extended by the following major 
assumption: For the compressible laminar boundary layer with heat trans- 

fer across the surface, the skin-friction and heat-transfer parameters 
I and r can be correlated only in terms of the parameters n and S^, 
derived from the exact solutions of reference 12. 

It is thus implied that the solutions of reference 12 adequately 
represent the general boundary layer, although they were derived for 
Falkner-Skan type flow. 

As in all first-order boundary-layer theories, the pressure distri- 
bution (and consequently the external velocity distribution) is assumed 
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to be known. Then the utility of the correlation may be stated as fol- 
lows: If n were known at a given point on the surface, 9 - (and hence 
e) can immediately be obtained from equation (22). If Z(n) is a known 
function for the specified wall temperature, the wall shear is immediately 
obtainable from equation (2l) . Similarly, if r(n) is known, the heat 
transfer can be found from equation (23) . 

If the postulate of correlation in terms of n and S be admitted, 
c/j equation (27) becomes 

oi 

- Ue ® j = H(n,S„) (26) 

vhere 

W(n,S^) = 2[n(H^^ + '2) + (29) 

^hls is the fundamental equation of the present method. Its solution, 
resulting in a determination of n(x), is the first stage in solving for 
the boundary-layer characteristics. Then the function J(n,S^) is used 

to determine the wall shear, and the function r(n,S„) is used to deter- 
mine the heat transfer. 


Evaluation of Correlation Parameters 


The quantities Z, n, and r defined in equations (2l) to (23), 
as evaluated from the solutions of reference 12, are listed in table II. 
An alternate parameter to r for the determination of the heat transfer 


is the Reynolds analogy parameter, defined as 


C^( 


w 


relating the heat 


Nu ^ 

transfer to the skin friction. Because this (arbitrarily chosen) param- 
eter is herein determined from solutions for a Prandtl number of 1.0, 

"fR^w') 

Nu A 


it will be denoted 


( 


'Pr=l 


as it appears in table II. It is related 


to 


by 


V Nu / 


23^2 


(30) 


Pr=l 


The parameters 
are plotted against 


\ Jft=1 


n and S__ are 

w 


n 


and I I as functions of 

/Pr=l 

bn figures 1 and 2, respectively. The solid 
portions of the curves represent the solutions of reference 12. The 
reversal of the curves for = 1.0 is associated with the velocity 
overshoot phenomenon discussed in reference 12. 
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When very strong favorable pressiire gradients are maintained for 
some distance along a surface, it may become necessary to select a cor- 
relation for a range of correlation numbers more negative than those 
resulting from the calculations based on reference 12. This may il- 
lustrated by comparison of the extreme case of a Falkner- an 
having an infinitely favorable pressure gradient (p = 2) with the case 
of a real flow having a fioced favorable gradient following a run of zero 
nressure gradient. In the extreme case it can be shown that the momen^.um 
thickness is always zero, but that the correlation number (proportiona- 
te -u^ 0^) is finite since u^^ is infinite. With the real flow, the 

momentum thickness in the region of favorable pressure gradient can be 
made as large as desired by increasing the length of run preceding the 
favorable gradient. Thus, the value of n may become as negative as 
desired. Under such circumstances extrapolated correlation flections 
would have to be used. The broken portions of the curves of figures , 

2, and 6 represent an extrapolation that is assumed to be qualitatively 
correct and quantitatively reasonable. 


METHOD OF SOLUTION 


The solution of equation (28) obviously depends on a knowledge of 
the term N(n,S ). This quantity was evaluated from equation ^,29j and 
the associated formulas for i, n, and The results are shown in 

figure 3. It is to be noted that the curve for = 0 is nearly a 

straight linej and for all negative values of the curves depart 

only slightly from this condition except in the range near separation 
wSe thf curves become double- valued. For > 0 (hot wall) the curve 

is essentially straight except in the region of strong favorable gradients 

The examination of these curves and of equation (28) suggests two 
methods of detemining the correlation number n. The first method i. 
that of solving equation (28) by numerical integration. This method .s 
assumed applicable to either variable or constant wall temperature. ..he 
necessary numerical procedure, however, is tedious, since it 
integration of a first-order nonlinear nonhomogeneous ordinary differ- 
SSE equation. Because a simpler method is available when the surface 
is isothermal, no numerical Integration procedure will be presented here 
in detail. However, some integration relations are stated in appendi 

The second method, applicable when the surface temperature is con- 
stant (or presumably, nearly constant), will be termed the linea 
method uses the nearly linear shape of the e^ves of 
against n for constant S„. It directly corresponds to the procedure 
of Thuaites for Incompressible flo« and to ttot of Kott and C rabtr ee for 
compressible flou over Insulated surfaces. The curve of H against 
for a given 


S is assumed represented by 
w 


N = A + Bn 


(31) 
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f 4-^^ equation (3l) is inserted in equation (28) a simple linear 
first-order equation results^ which has for its solution 


n = - AU, 


-B 


n 


/ 


-X ^0 


B-l 

Ue dX 


(32) 


If equation (32) is transformed to physical quantities by usine 
Stevartson s transformation^ there results for two-dimensional flow 


n 


^0 

t^ 


(%r f 


|L /t„\K 




M 


B-l 


^0 


>(f) 


(33) 


-1 

where K — -l) ■’ ^ fixed leng)th^ and the dimensionless pres- 

sure gradient P' is given by 


dx 


L 


dx 


Pe^e^ rPgM^ 


(34) 


member of actuation {33) has been arranged in a form convenient 
for later use. 


The determination of the coefficients A and B is as follows: 

If the straight line (31 ) is chosen to pass through the correct value 
of N at zero pressure gradient (n = O), then A « 0.44 independent 
of S . In this case^ only B is affected by the presence of heat 
trans er . Figure 4 shows the values of B(S^) for the following choices 
of matching conditions: 


(l) The line (31 ) goes through the point corresponding to N = 0 
^Ue = O) in order to match the conditions at a stagnation point. Thus, 


sp 


* / CLL» ex 

for two-dimensional flow (n = o)^ B = - 0.44/n 

(2) The line (31 ) coincides with the tangent to the u(n) lines at 
n = 0 (small pressure gradients). 

(3) The line (3l) is selected to give good over-all agreement for 
unfavorable pressure gradients . 


If a better match with the curves of figure 3 is desired in calcu- 
lating n for certain ranges of pressure gradient, a tangent line to 
the curve of K against n may be chosen at a desired value of n. 
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For instance, in considering tlie flow in the vicinity of a two-dimensional 
stagnation point with = -0.8, the tangent line through K = 0 has 

the coefficients A = 0.372, B = 2.53. _ The value of B is ^uite sensi- 
tive to the matching assumption, especially in the region ^ 0 ^iiS- 

4). However, the final value of n(x) is somewhat imsensitive to ihe 
value of B, since the terms involving B in equation (,33) appear both 
inside and outside the integral in a compensating manner. The acc^acy 
of the method decreases, of course, in regions where the plots of iig- 
lire 3 have large curvature . 

The calculation procediire is as follows: Values of A and B are 

chosen for use in equation (33) either from figure 4 or from tangent- 
line considerations. The integration is then perfoimed by us^g a suit- 
able integration rule and a proper step size. It is reco^ended that 
the step size chosen be as small as practicable in order to obtain re- 
sults which are reasonably smooth. In some cases (e.g., near e 
boundary-layer origin) it may be advisable to perfoim the _ integration 
by obtaining Taylor's series expansions of the integrand in the van 
(x/l). Then the integration can be carried out in closed form, corre- 

sponding erfcdively lo zero step size* 

There are two possible starting conditions in a boundary-layer cal- 
culation: (l) that of a sharp edge, that is, 0 = 0, n - 0; or (2) he 
sS^ation poLt, where Ue = u, = 0. In using the linear method, the 
starting conditions are automatically satisfied when the chosen line 
(31) goes through the starting point. Thus, if matching condition (1) 

Is used, both possible starting conditions can be satisfied, since the 
corresponding line (3l) passes through the curve from the exact so u- 
tions at both n = 0 (N = 0.44) and, for the two-dimensional flow at 
N = 0 (n = -0.44 /b). Values of n for stagnation-point flow taken 
from fig^e 3 are shown in figure 5. 

The corresponding relations and procedure for axially sy^etric 
flow based on Mangier 's transformation (ref. 26) are presented in ap- 
pendix C . 

It is sometimes helpful to have an analytical expression for the 
initial variation of correlation number as a check on the numerical cay 
SSaiSns Se Initial variation ot n witn x the various start- 

ing conditions, as represented hy the derivative (s)^' 
in appendix: D. 


3321 



3321 


MCA TN 3326 


15 


BOUNDARY-IAYER CHARACTERISTICS 


Once the correlation nninber n is detertained as a function of 

it is possible to obtain I and — j from figures 1 and 2. 

Then, the local skin-friction coefficient and heat transfer are easily 
calculated from the following relations, which apply to both two- 
dimensional and axially symmetric flows: 

XI 


From the definitions C^. = — 


w 


UeX 


-, Re„ = , Nu = -r 

>■’ w V ' t 


(IX 


2 


w 


- t 
a^w w 


and 


from equations (2l) and (22)^ it follows that 



(35) 


or that (the equivalent form which may be more convenient) 


CfV^e^ 



(36) 


Once C^ is determined, the heat transfer may be calculated from 
curves of the Reynolds analogy parameter against correlation number of 
figure 2, by using the relation 


Nu 




(37) 


In utilizing equations (36) and (37), it is useful to have the 
initial values of the parameters. These values are listed in table I. 
It may be noted that, at a stagnation point, equation (36) reduces to 

Cf-^Re^ = (38) 

The calculations thus far have been for Pr = 1.0. The effect of 
Prandtl number on skin friction is small and is therefore usually neg- 
lected. It can be seen from table I that, for stagnation-point flow, 
the maximum difference in the quantity Cp-y/Re^ between solutions 
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for Pr = 1.0 and Pr = 0.7 is about 7 percent. With regard to heat 
transfer, Tifford and Chu (ref. 27) have fo\md, from solutions with con- 
stant fluid properties, that the effect of Prandtl number on heat trans- 

/ \ / 

fer can be accounted for by multiplying ( y ■ ) by (Pr) . Values 


\ V W/pp_p 

of a suggested in reference 27 are as follows: For small pressure 

gradients, a = l/3; for large adverse pressure gradients, a = l/4; and 
for extreme favorable gradients, oc = l/2. Squire (ref. 28) has indicated 
that a = 0.4 is adequate for stagnation-point flows. Recently obtained 
solutions (as yet unrepcrted) of equations (16) of reference 12 for 
P = 1, Pr = 0.7, and Mg ® show that this type of correction may be 
adequate for all compressible boundary-layer calculations. 


It should be noted that, in the definition of Nusselt number, the 
temperature difference in the denominator was assumed to be (tg^^^ - t^) . 
Since for Pr = 1 the recovery temperature is tg, the present calcu- 
lations (based on those of ref. 12) can give no indication of the adia- 
batic wall temperature for Pr 1. For a first approximation, it may 
be reasonable to calculate t„ „ by using a temperature recovery factor 

of (Pr)^/^. This is the well-known expression for recovery factor for 
the case of high-speed flow with zero press\nre gradient. The accuracy 
of its application to flows with large pressure gradients is not yet 
established. 


The physical momentum thickness is determined from 



(39) 


The displacement thickness 5* may be calculated by using the fol- 
lowing simple expression for the ratio of displacement thickness to m.o- 
mentum thickness: 



(40) 


In reference 22, this expression was suggested for use in calculating 
flows over insulated surfaces with Pr = 1. For noninsulated surfaces, 
the dependence of upon wall temperature is presented in figure 6. 

With large amounts of cooling in favorable pressure gradient flows, it 
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factors (here corresponding to negative dis- 
placement thicknesses) result. This occurs because tSe surffce cLlSg 
produces an increase in density near the wall, so that there is more 

^tLmrfloi."^^^ boundary layer than in the 


CO 

H 


EXAMPLES 


An important test of the method 
final results for practical problems 
or vith experimental results. 


developed is the comparison of the 
the findings of other theories 


I. Supersonic Surfaces 


The linear method for determining the correlation number n is 
applied to the calculation of skin friction and heat transfer for the 
two supersonic surfaces at Mach number 3.0 calculated in reference 6. 
These surfaces are shown in sketch (b) . 



(b) 

A comparison is made, in the following table, between the results 
obtained by using the linear method and those obtained by using Low's 
perturbation method with Pr = 0.72 (ref. 6). The comparison is made 
for a hot wall, an insulated wall, and a cold wall at x/L = 1 (see 
sketch (b)). A value of a of 1/3 was used in these calculations. 
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Gradient 

Heat -transfer 


(Hu/ 

coiidxtioii 

Linear 
method, 
Pr = 1 

Low (ref. 6), 
Pr = 0.72 

Linear 
method^ 
corrected 
for Pr 

Low 

(ref . 6 ) , 
Pr = 0.72 

Favorable 

Cold wall 

0.680 

0.679 

0.270 

0.271 


(S^=-0.9) 


.876 




Insulated 

.877 




wall 






Hot wall 

1.031 

1.048 

.284 

.318 


(S^O.61) 





Adverse 

Cold wall 

0.645 

0.661 

0.320 

0.307 


(S^=-0.9) 


.475 




Insulated 

.406 




wall 






Hot wall 

.256 

.386 

.274 

.277 


(S„=0.43) 






A comparison of values indicates agreement of .^kin friction ^if^^ 

2 nercent in the case of a favorable pressure gradxent. For the adve.s 
prLsure gradient cases, reasonable agreement is obtained for the cold 
Sail, allhough for the insulated and hot walls, th.^ 

obtained. This difference is due essentially to the f^^f ' 

case of an adverse pressure gradient, the solutions of C ' 

which resemble the series -type solution fj^f 25) 

from Falkner-Skan type solutions such as that ’ble is 

(p p sketch (a)). Some of the difference in the preceding table is 

effect; however, for the Seated s^fa« an adveree 

nressure gradient, the flow is closer to separation than appears per 
.issible for a theory based be 

:;“LbL! SoH^riLnt is also obtained for heat transfer except 
for the case of the heated surface with favorable pressure gradient. 

Some of that difference might be a Prandtl number eflect. 

II. Axially Syimnetric Convergent -Divergent Rocket Nozzle 

The second example, that of a rocket nozzle, is one 
laree nressure gradients and*heat transfer. The nozzle chosen is illus- 
trated in figurf 7. It has a 25° half -angle convergent section ^d a 

Lir&Tr0CL XII XX^UJ-c: i T 4 - f Vi ^Tn'KoT' ciX.ft crma.'t', i DTI 


15 " 


hrarl in fi (TUre 7. fU UaS a ZO , . „ 

“^“angfTdlvergent section. The conl>ustlon-chamber stagnation 
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pressure is assumed to be 500 pounds per square inch absolute, the stag- 
nation temperature is taken as 4000° R, and the Prandtl number is assumed 
to be 0.78. The nozzle wall is assumed cooled to a uniform temperature 
of 800° R, which corresponds to S^. = -0.8. For the assumed 3-lnch 
throat diameter, the rocket has a nominal thrust of 5550 pounds for 
Y = 1.3. Local static conditions along the nozzle wall were obtained 
using one -dimensional flow relations. 

The calculation was performed by the linear method with A = 0.372 
and B = 2.53. In order to eliminate the effect of step size in the in- 
itial portion of the integration, a series expansion of the integrand was 
used for 0 < (x/l) < 0.5. For (x/l) > 0.5, the step size taken was 0.1. 

The resulting variation of n in the nozzle is also shown in figure 7 . 

It is seen that, in a portion of the nozzle including the throat 
(l.2 < (x/d) < 3.5), values of n are obtained which are more negative 
than those obtainable from the exact solutions. Thus, in order to cal- 
culate skin friction, heat transfer, and displacement thickness from 
1.2 <(x/D) < 3.5, it is necessary to use the extrapolated portions of 

the curves of figures 1 and 2 for = -0.8. Wo extrapolation is 

needed to obtain momentum thickness since the momentum thickness is re- 
lated to n through equation (39). 

The calculated local heat -transfer rates as well as displacement 
and momentum thickness are shown in figure 8. It is seen that large 
rates of cooling are required in the neighborhood of the stagnation 
point and the nozzle throat. If the cooling were to become insufficient, 
these seem to be the most likely locations of failure. The required 
local cooling rate to maintain constant wall temperature decreases 
sharply beyond the throat of the nozzle. 

In the absence of more appropriate information, it has been cus- 
tomary in recent years to use flat— plate heat— transfer relations in 
estimating heat transfer in a nozzle. The use of such a relation in 
the current problem (Wu/ -^Re^ = 0.305 for Pr = 0.78) yields values 
indicated by the dashed line in figure 8. It is seen that for this 
problem the flat -plate relation seriously underestimates the amount of 
cooling required over a large part of the nozzle . This illustrates the 
importance of considering the effects of pressure gradient on heat 
transfer . 

The momentum thickness is seen to reach a minimum value at the 
nozzle throat. The displacement thickness is a small positive quantity 
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at the stagnation point but is negative for most of the convergent sec- 
tion as well as in the vicinity of the throat.^ 

The use of different values of A and B in performing the mmer- 
ical integration would have approximately the following effect: With 

the values A = 0,44 and B =* 3.0 (from fig. 4)^ the momentum thick- 
ness would be about 10 percent smaller in the vicinity of the throat 
than the values in figure 7 and would be within 5 percent of the pre- 
sented values over the rest of the nozzle. With A = 0.335 and 
B = 2.34^ the momentum thickness at the throat would be about 6 percent 
larger than the value presented. The effects of varying A and B on 
skin friction and heat transfer would be less than 3 percent at the 
throat . 


CONCLUDING REHyiARKS 

The application of Stewartson's transformation to the con^jressible 
laminar bo\mdary- layer equations with heat transfer yielded a siir5)le 
first-order system of ordinary differential eqiiations^ the first of 
which is very similar to the Karman momentum integral. Dimensionless 
shear and heat -transfer parameters were defined. The assumption of 
correlation of these parameters in terms of a momentum parameter resulted 
in a complete system of relations for calculating skin friction and heat 
transfer. Knowledge of velocity or temperatiire profiles is not neces- 
sary in using this calciilation method. Procedures for the calculation 
of the longitudinal distribution of correlation number are presented^ 
which include as special cases the method of Thwaites and that of Rott 
and Crabtree. The dimensionless parameters introduced herein were 
evaluated from the exact solutions of reference 12. 

Calculations of an example involving small pressure gradients show 
that the method is reliable when compared with the perturbation method 
of reference 6 over the same range of Mach number^ pressure gradient^ 
and heat transfer. 

The method is also applied to the calculation of heat transfer and 
displacement thickness in a highly cooled^ convergent -divergent^ axially 
symmetric rocket nozzle. The results of this calculation show that high 
rates of heat transfer are obtained at the initial stagnation point and 


3 

This unusiml result produces the interesting possibility that^ 
for a rocket nozzle with cooled walls and viscid flow, a mass discharge 
coefficient based on throat area, generally assumed to be less than 
unity because of boundary-layer "blockage” at the throat, may actually 
exceed 1. A distinction exists between this phenomenon and that of 
negative momentum thickness (refs. 1, 10, and 12) associated with 
"velocity overshoot . " 
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at the throat of the nozzle. Also indicated are negative displacement 
thicknesses in the convergent portion of the nozzlej these occur be- 
cause of the high density within the lower portions of the cooled 
boundary layer. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio, February 1, 1955 
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APPENDIX A 
SIMBOLS 


A 

a 

B 


The following s^hols are used in this report: 
constant from N = A + Bn 
sonic velocity 
constant from N = A -i- Bn 


local skin-friction coefficient^ 


2T 


w 


P U 
w e 


D 

E 


specific heat at constant pressure 
nozzle-throat diameter 


enthalpy-flux thickness ^ E 



S 


H 

H 

h 

K 

k 

k 


inc 


su 

L 

2 


stagnation-enthalpy-defect thickness^ S 



form factor^ H = 

physical form factor for - 0^ ^inc 


enthalpy 

37 • 1 
2(7 - 1) 

thermal conductivity 




Sutherland *s constant 


arbitrary length 


dimensionless shear parameter, 2 



dY 
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M 

in 

N 


Nu 


n 




Pr 

P 

q. 

E 


Mach numher 


exponent from Falkner-Skan external velocity distribution 

Uq = cx™ 


momentum parameter 


, N = 2 r 


^^%nc + 2) + 2 


Wusselt number, Nu = 


\^)m 


- t 


■] 


w 


correlation number^ n 



dimensionlesB pressure gradient^ P^ 



2 

© 


Prandtl number, Pr = — 2- 
^ k 

static pressure 

local heat-transfer rate^ Btu/(sq. ft) (sec) 
radius of axially symmetric body 


Ee 

w 


P U X 

Eeynolds number, Ee = — ^ ^ 


r 

S 

t 

t 

a^w 


heat -transfer parameter^ 



h 

enthalpy function, S = r-2- - i 
’ h 

0 


static temperature 
adiabatic ■wall temperature 
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transf onusd longit.udinal velocltiy ^ U u. ^ >|*Y 

0 

longitudinal velocity component 
transformed normal velocity, V — ~ 
normal velocity component 

transformed coordinate along surface, X = 


r 


X ^ dx 

Pq ^0 


X 

Y 

7 

a 

P 

7 

5* 

9 

\ 

V 

P 

T 


coordinate along surface 
transformed normal coordinate^ Y 

normal coordinate 

exponent of Prandtl number in Eeynolds analogy parameter 

2 m 

pressuro gradient parameter, ^ m + 1 
ratio of specific heats 
displacement thickness 
momentum thickness 



dynamic viscosity 
kinematic viscosity, v = p/p 
mass density 
shear stress, t = P 3 — 

07 

stream function, eq. ( 6 ) 
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03 

03 

ro 


Subscripts: 

e local flow outside boundary layer (external) 

i associated transformed q.uantity 

B local stagnation value 

sp stagnation point 

w wall or surface value 

0 free -stream stagnation value 

1 initial value 
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APPETTOIX B 


MJMEEICAL INTEGEATION METHOD 

Method. - The most direct method of solving equation (28) is 

numerical integration, using the calculated curves of N(n,S^) for deter 
mination of the right member. This method is assumed applicable to 
either variable or constant S^. 

An integration procedure may be simply Indicated by direct inte- 
gration of equation (28). The resultant Integral equation can be 
written : 


For two -dimensional flow, 


n 





(Bl) 


For an axially symmetric closed body, through Mangier 's transformation 
(ref. 26), 


n = 




(B2) 


Since the integrands contain N(n), which is unknown at this point, 
no simple evaluation is possible. In fact, these equations are actually 
only a condensed notation of the procedure to be followed. The Inte- 
gration must be carried out piecemeal, alternating with determination 
of the left member of the equation and iterating for accuracy. 


The necessity of working with the transfonned coordinates can be 
eliminated by considering the Stewartson transformation from to u^. 

For example, in physical coordinates equation (B2) becomes 
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where L is any fixed length, and the dimensionless pressure gradient 
P ' is given by 


<iPp * dp 
^ 

p, ^ ^ ^ ^ (B3) 

%e4 'rPeM^ 

If the Isentroplc relation p/p^ = constant is used in equations 
(31) and (32) and if the Stewart son transformation is applied^ there 
results : 


For two-dimensional flow, 





For an axially symmetric flow (closed body), 


X 



Initial values . - When the numerical integration method is used, 

certain considerations are necessary in order to start the solution 

properly. There are two possible starting conditions: (l) sharp edge 

or pointed body, where 0=0 and n = 0, and (2) stagnation pointy 

where U = u^ = 0 . 
e e 

In the case of a boundary layer starting from a stagnation point, 
the initial value n^ of n is determined from the condition = 0, 

Uo = constant in equation (28). For two-dimensional stagnation-point 
© 

X 

flow, the Hartree pressure gradient parameter P is equal to 1.0. 
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Since^ for the Falkner-Skan type flow considered, N 


3 


n, 


It is 


seen that N = 0 at a two-dimensional stagnation point. This fixes 
n at the values shown in figure 5, which were obtained from figure 3. 
For axially symmetric stagnation-point flow over a closed body, p = l/2 
(ref. 29), so that = -2n^. The values of n^ for axially symmetric 

stagnation-point flow over a closed body as obtained from figure 3 are 
also shown in figure 5. For the stagnation-point flow over the blunt 
lip of an open axially symmetric body, n^ can be shown to have the 
two-dimensional value. 
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APPENDIX C 


LINEAR METHOD FOR AXIALLY SYMMETRIC FLOW 


For axially symmetric flows, the following eq.uation, which Is 
equivalent to equation (33), Is obtained by application of the trans- 
formation of Mangier (ref. 26): 


K 1-B f /+ 

_n_ ^ /M R^ 

I 2 / Vt^l „ 1-B 

■ -0 \ 0>' R J \ OJ % 


>(!) 


(Cl) 


t" 

e 


where E = E(x) Is the radius of the hody at station x, K = 


P' = 


dx 


2 

7P M 

e 


In evaluating the coefficients A and B of the straight line 

N = A + Bn (31) 


A may he chosen as 0,44 so that the line (31) passes through the cor- 
rect value of N at n = 0. The choice of B may he made so that the 
line (31) goes through the point where N = -2n in order to match the 

( 0 44 
— ^ + 2 

Pep ^ 

For achieving better agreement with the curves of figure 3 in certain 
ranges of pressure gradient, a tangent line to the curve of N against n 
may he chosen as was indicated for two-dimensional flow. 

The three possible starting conditions are: (l) For a pointed 

hody, 8 = 0, n = 0. (2) For a stagnation point on a closed body, 

^e ” ^e “ ^ that N = -2n, Values of n^p for this axially 

B^ymmetric stagnation point are shown in figure 5 and indicated in 
table II for p = l/2. (3) For a stagnation point on the blunt lip of an 
open axially symmetric body, it can be shown using the axially symmetric 
form of equation (28) that N = 0 so that n is that for two- 
dimensional flow. 
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APPENDIX D 


INITIAL VARIATION OF CORRELATION NUMBER 

It is sometiznes helpful to have an analytical expression for the 
initial variation of correlation number as a check on the numerical 

calculations. The following expressions for 
from equation (28). 


are determined 


x=0 


Sharp edge : 


Two-Dimensional Flow 


= - 0.44 f^o] 

^^<=0 % \~eJ 

stagnation point (blunt body) : 


x=0 1 + 


^sp r«xx\ 

ifrw 

' 'sp sp 


(1)1) 


(D2) 


Axially Symmetric Flow 

For axially symmetric flow the initial derivatives must be evaluated 
from the following equation: 


^ A. 

r2 


= K(n,S^) 


(D3) 


Closed Body 

Pointed nose : 


= *. 0.147 po^ 

V'^^>'x=0 VW 


(D4) 


3321 



3321 


MCA TN 3326 


31 


Stagnation point (blunt nose) : 



Open Body 

Sharp lip: 



Stagnation point (blunt lip): 
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TABLE I. - INITIAL VALUES OF PARA.METERS 


(a) Stagnation-point flow 


Sw 

Cf 


Nu/V^ 



H , 
sp 


Pr = 1 



Pr = 0.7 
(b) 


Pr = 1 

Pr = 1 

Pr = 1 


Two-dimensional (p = l) 


-1,0 

'^1.30 

1 .21 

^0.506 

0.458 

0.436 

^-0.170 

^0.400 

^-0.425 

-.8 


1.49 

^ .522 

.452 

.451 

c .012 

^ .384 

.031 

-.4 

^2.04 

2.00 

.546 

.474 

.475 

'= .345 

^ .338 

c 1.021 

0 

2.46 

2,46 

.570 

.495 

.495 

.648 

.292 

2.218 

1.0 

3.47 

3.54 

.615 

,553 

.537 

1.386 

.177 

7.850 





Axially symmetric 

(P = 1/2) 



-1.0 

1.64 

1,56 

0.700 

0.607 

0.607 

-0.0771 ] 

0.300 ^ 

-0.257 

-.8 

1.85 

1.79 

.712 

.617 

.621 

.0576 

.289 

.199 

-.4 

2.25 

2.22 

.739 

.639 

.643 

.318 

.269 

1.165 

0 

2.62 

2.62 

.763 

.662 

.662 

.569 

.248 

2.298 

1,0 

5.49 ' 

3.55 

.809 1 

.701 

.708 

1.165 

.194 

6.012 


^hese values 
'These values 
'^'interpolated 


are obtained when eqs. (16) of ref. 12 are solved for Pr = 0.7, 
are obtained multiplying the results for Pr = 1 by 
values from solutions of ref. 12. 


Me = 0. 


(b) Sharp edge or pointed body 


CfV^ 

Nu/-^/Re^ 

Pr = 1 

Pr = 0.7 
(a) 

Tw o - d i men s i onal 

0.664 

0.332 

0.295 

Axially symmetric 

1.150 

0.575 

0.510 


^These values are obtained 
by multiplying the 
results for Pr = 1 by 


(0.7) 


1/3 _ 
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TABLE II. - SUMMARY OF HEAT-TRANSFER AND WALL-SHEAR PARAMETERS 



6 

n 

1 

N 

r 

/CfRe\ 

VNu 

%nc 

-1.0 

-0.326 

0.1335 

0 

1.0845 

0.0212 

0 

2.063 


-.3657 

.1579 

.3929 

1.1804 

.0307 

.3381 

1.530 


-.3884 

.1591 

.0896 

1.1382 

.0359 

.7939 

1.013 


-.360 

.1257 

.1446 

.9504 

.0297 

1.224 

.630 


-.30 

.0907 

.1749 

.7858 

.0212 

1.493 

.404 


-.14 

.0343 

.2063 

.5590 

.00774 

1.830 

.134 


0 

0 

.220 

.440 

0 

2.000 

0 


.50 

-.0897 

.2459 

.1793 

-.0188 

2.347 

-.257 


2.00 

-.2938 

.2829 

-.2933 

-.0586 

2.837 

-.538 

-0.8 

-0.3088 

0.1215 

0 

1.0305 

0.0172 

0 

2.240 


-.325 

.1304 

.0312 

1.0606 

.0210 

.3100 

1.828 


-.3285 

.1298 

.0436 

1.0499 

.0216 

.4194 

1.708 


-.3285 

.1260 

.0681 

1.0185 

.0220 

.6245 

1.501 


-.325 

.1212 

.0827 

.9885 

.0216 

.7438 

1.396 


-.30 

.1017 

.1214 

.882 

.0187 

1.058 

1.138 


-.14 

.0355 

.1935 

.5781 

.00642 

1.712 

.692 


0 

0 

.220 

.44 

0 

2.000 

.518 


.50 

-.0837 

.2678 

.1676 

-.0138 

2.599 

.199 


1.50 

-.2008 

.3179 

-.1332 

-.0313 

3.263 

-.083 


2.00 

-.2522 

.3366 

-.2517 

-.0388 

3.502 
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Figure 3. - Correlation of dimensionless momentum equation. 
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Enthalpy function at wall^, 

Figure 4. - Variation of B for use in linear method of 
determining correlation number. 









Figure 6. - Form-factor correlation. 
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Figure 7. 


- Variation of correlation number in rocket nozzle. 
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Figure 6, - Results of rocket-nozzle calculations. 
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